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In these lecture notes I re- visit the problem of the deep number counts as a viable 
method to (1) study the evolution of the galaxy population in the past and (2) determine 
the value of the deceleration parameter go- After a brief description of the basic concepts 
needed to model the counts, we show that simple luminosity evolution models which 
incorporate dust are able to provide reasonable fits to both the deep counts and the 
redshift distribution of faint field galaxies if go = 0.05. In order to fit the counts in a 
closed (qo = 0.5) Universe, we include a population of fading dwarfs, more numerous in 
the past, based on current ideas to understand the nature of dwarf galaxies. 



INTRODUCTION 



Galaxy number counts have been widely-used in Cosmology 
in an attempt to estimate the value of the deceleration pa- 
rameter go (throughout this paper I will assume A = 0). 
^— ^ Whereas most observational approaches are based on the 
direct estimation of the density of matter in the Universe 
Z~_'i (es. superclusters mass/volume), number counts deal with 
^3 the pure geometric effect of varying go . 

The idea is fairly simple. Let us imagine that we take a 
very deep wide-field image at any place on the sky, counting 
afterwards the number of galaxies per unit area and magni- 
^ 1 tude bin (7V(m)). The value of go can be inferred by means 
O ipf ihp comparison of the observed N (m) with the predictions 
^from different go models. In the most basic models, usually 
» I called no-evolution models, two assumptions are made: (1) 
the luminosity of the galaxies remains constant on time, (2) 
the number of galaxies per unite volume is conserved after 
they were formed. 

Since the first papers on number counts (Brown & Tins- 
ley 1974; Bruzual & Kron 1980; Shanks et al. 1984; King & 
Ellis 1985; Koo 1986) it has been well known that these 
simple no-evolution models are unable to fit the counts, no 
matter what value of go is used. There is an excess of galaxies 
at faint magnitudes (B fainter than ~ 22 — 23 depending on 
the normalization), which is much larger when the counts are 
made using blue photometric passbands. Thus, this problem 
is usually referred to as the excess of faint blue galaxies. 

In order to reconcile observations with model predic- 
tions it is neccesary to account for evolution in the galaxy 
population. The luminosity of the galaxies is expected to 
vary with time, as stars are evolving. In particular we ex- 
pect that galaxies were brighter at earlier times, closer to 
the epoch at which they were formed. After the first episodes 
of star formation the birth rate of stars is expected to de- 
crease (either because the interstellar gas is exhausted or 
because it is blown out of the galaxy), like it seems to be 
the case in early type galaxies, or to remain approximately 
constant on time, as we observe in late type galaxies. If 



galaxies were brighter in the past then for a given apparent 
magnitude they are observable at greater redshifts, and so 
for the same angular area we survey a larger volume (Tins- 
ley 1980; Guiderdoni & Rocca-Volmerange 1990; Metcalfe 
et al. 1991). 

It has also been suggested (Guiderdoni & Rocca- 
Volmerange 1991; Broadhurst et al. 1992) that not only the 
luminosity but also the number density of galaxies could 
vary with time. This variation is, for example, expected in 
hierarchical models of galaxy formation where galaxy sub- 
units merge to build-up the present day population. In these 
models galaxies are assumed to be more numerous in the 
past, and so they predict larger numbers of galaxies to be 
seen at faint magnitudes. 

In this paper we re-visit the modelling of the number 
counts by means of luminosity evolution models which in- 
corporate dust. In Section 2 we briefly summaryze the basic 
concepts to model the predictions in the no-evolution case. 
The effects of the luminosity evolution (and the dust) in the 
model predictions are analysed in Section 3. Section 4 is de- 
voted to the study of the redshift distributions of faint field 
galaxies, as an independent test to the number counts mod- 
els. In Section 5 we show that in order to fit the counts in a 
go =0.5 Universe in the framework of luminosity evolution 
models, we must allow for luminosity-dependent luminosity 
evolution. We also show that the introduction of a popula- 
tion of fading dwarfs allows to fit the counts. Finally, the 
conclusions are presented in Section 6. 



2 MODELLING THE NUMBER COUNTS 

The number of galaxies with apparent magnitude [m,m -\- 
dm\ located in the redshift shell [z, z + dz], which contributes 
to the counts per steradian can be computed as follows: 



, , dV 

d^A(m,z) = $(M)(1 + zf — dmdz 

dz 



(1) 
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Table 1. Luminosity Functions 



Galaxy type 


</.♦(%) 


a 


M* 


Early 


57 


-0.7 


-19.6 


Intermediate 


26 


-1.1 


-19.9 


Late 


17 


-1.5 


-20.0 



where $(M) is the luminosity function (LF; it can be esti- 
mated using nearby redshift surveys) and M is the absolute 
magnitude of a galaxy observed to have an apparent magni- 
tude m when is located at a redshift z. dV/dz is the volume 
element per steradian of the redshift shell [z,z + dz], which 
in the particular case of Friedman models depends on the 
values of go and Ho (see eg. Weinberg 1972). 

The relation between the absolute magnitude of a 
galaxy with its apparent magnitude at any z is given by: 

m = 5\og{DL} + M + K{z) (2) 

where Dl is the luminosity distance, which again in Fried- 
man models depends on the values of go and Ho. K{z) is 
the /('-correction, which accounts for the red-shifting of the 
spectral light due to the Doppler effect. For any given pho- 
tometric passband we observe different wavelength ranges 
(in the rest-frame of the galaxy) depending on the redshift 
at which galaxies are located, and so we have to do the ap- 
propiate correction because the spectra of galaxies are not 
flat. For the no-evolution models, the /('-corrections are com- 
puted using spectra of nearby galaxies. 

The number of galaxies per steradian and magnitude 
bin can be obtained by integrating d^ A from 2 = up to 
the redshift at which galaxies were formed Zfor, 

N{m)dm= I d^A{m,z) (3) 

J z = 

In what follows we will use the Schechter parametriza- 
tion of the LF, which has the form: 

$(A/) =0.92$'exp[-0.92(a-M)(A/- A/') 

- exp[-0.92(A/ - M*)]]dM (4) 

where a is the slope of the function at its faint end and M* is 
a characteristic magnitude at which the knee of the function 
occurs. 4>* is the normalization, i.e. the number density of 
galaxies. Notice that both 4>* and M* depend on Ho- The 
dependence of 4>* and dV/dz on Ho cancels out in equation 
(1). On the other hand, the dependence of Dl on Ho cancels 
with that of A/*. Therefore the number counts only depends 
on the value of go (and of course A, in the case it had a value 
different than zero). 

In Figure 1 we show the number counts in the B-band 
together with predictions from no-evolution models. The 
data is a compilation from the literature: Kron 1978; Koo 
1986; Jarvis & Tyson 1981; Couch & Newell 1984; Infante et 
al. 1986; Tyson 1988; Jones et al. 1990; Lilly et al. 1991; Met- 
calfe et al. 1991, 1995. To compute the predictions we have 
assumed three different galaxy types, each one with different 
/('-corrections and LF. The LFs (see Table 1) have been ob- 
tained from the Durham-redshift survey (Shanks et al. 1984; 
see also Metcalfe et al. 1991). As we already mentioned in 
the introduction, the no-evolution models are unable to fit 
the counts, even for low-go values for which the volume at 
high- 2 is larger. 
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Figure 1. The number counts in the B band. The data has been 
taken from the literature (see text for references). Also we plot 
predictions from no-evolution models. 

3 LUMINOSITY EVOLUTION MODELS 

The luminosity of the galaxies is expected to evolve with 
time. A simple method to include luminosity evolution in 
the models is by adding an extra-term in the right-hand 
side of equation (2) as follows: 

m{z) = 51og(/)i,) + M(z = 0) K(z) + e(z) (5) 

where e(z) is the evolution correction, which measures the 
variation of the luminosity (in units of magnitudes) from 
redshift 2 to 2 = 0. In this equation we now relate the ap- 
parent magnitude of a galaxy located at a redshift z with 
the absolute magnitude that the same galaxy will have at 
2 = (this is because in equation (1) we will make use of 
the LF estimated using nearby surveys, i.e. at 2 = 0). 

The /('- and e- corrections (/(' -|- e in what follows) are 
computed by means of spectrophotometric models of stel- 
lar populations (eg. Tinsley 1972; Bruzual 1981; Arimoto & 
Yoshii 1986, 1987; Bruzual & Chariot 1993). These models 
are based on stellar libraries which provide information on 
the stellar spectra (or colours) at any position in the HR dia- 
gram, i.e. at any stage of the stellar evolution. We can model 
the spectral evolution of a galaxy by assuming a redshift of 
formation (z/or), an Initial Mass Function (IMF) for the 
stars and a Star Formation Rate (SFR). To check whether 
the chosen model parameters (z/or, IMF and SFR) are ap- 
propiate one compares the synthetic spectrum at 2 = with 
observed spectra of nearby galaxies. Finally, once we have 
the spectra of the model galaxy at any redshift, it is straigh- 
forward to compute its magnitude (at any z) just simply 
by convolving the (red-shifted) synthetic spectrum with the 
photometric passband. The difference between the magni- 
tude at z and at 2 = will give the /(' -|- e corrections. 




Figure 2. K-correction plus evolution for three types of galaxies. 
The age of the galaxies at ^ = is 16 Gyr., which corresponds to 
a redshift of formation of ^ = 6 for go = 0.05 and Hq = 50 km 
s-1 Mpc-1. 



We notice that when the luminosity evolution is in- 
cluded to model the counts, the predictions depend (~ 
weakly) on the value of Ho . The reason is found in the re- 
lation between age and redshift (which depends on Ho) to 
trace the galaxy evolution. 

In Figure 2 we have plotted the K + e corrections for the 
B-band for three model-galaxy types: early (E/SO), interme- 
diate (Sa/Sb/Sc) and late (Sd/Irr). The age of the galaxies is 
16 Gyr (zf or = 6 for go = 0.05 and Ho = 50 km s"^ Mpc"^) 
in all cases, and we have used the Salpeter IMF. The SFR 
is assumed to decay exponentially, with different e-folding 
times: 0.5, 7 and 900 Gyr for the early, intermediate and 
late types respectively. To compute the K + e corrections 
we have used the spectrophotometric models of Bruzual & 
Chariot (1993). As can be seen in the Figure, the K + e cor- 
rections for late-type galaxies are very small, because their 
SFR is nearly constant but also because their spectra are 
quite flat. On the contrary, for early type galaxies the K + e 
corrections are very large. The observed luminosity in the 
B-band for 2 > is smaller than for 2 = 0, up to a redshift 
of 2 ~ 2. The reason is that the UV-light enters into the B- 
band, and early-type galaxies are not strong emitter in the 
UV. However for redshifts z > 2 the galaxies become more 
and more luminous, as we move closer to the epoch were the 
bulk of stars was formed. 

In Figure 3 we plot the observed number counts in the 
B and K bands together with predictions from go = 0.05 
and go = 0.5 luminosity evolution models. The data in the 
K-band has been taken from Mobasher et al. 1986; Glaze- 
brook et al. 1993; Jenkins & Reid 1991 and Gardner et al. 
1993. To compute the models we have assumed three dif- 



Figure 3. The number counts in the B and K bands, and pre- 
dictions from luminosity evolution models with dust. 

ferent galaxy types as in Figure 2, and the LFs shown in 
Table 1. For the go = 0.5 model we assume 12.7 Gyr, which 
corresponds to Zfor = 10, rather than 16 Gyr for the age of 
the galaxies (the reason is that in a go = 0.5 and Ho = 50 
model the age of the Universe is only 13.1 Gyr). The models 
have been normalized to fit the number counts in the B-band 
at 5 = 17.5, which corresponds to an overall normalization 
of ({)* = 1.4 X 10"^ h"^ Mpc"^ The counts in the K-band 
were derived using the luminosity function in the B-band by 
means of the B — K colours computed for each model galaxy. 

To model the number counts we have also taken into 
account the presence of dust in galaxies. The effect of dust, 
neglected in many works on number counts, is proving to be 
very important (Wang 1991; Campos & Shanks 1995). The 
dust absorbes the light coming from the stars, re-radiating 
it in the IR. In general the effect of dust is much more severe 
as the wavelength shortens. This fact is reflected in the ex- 
tinction law, which gives the dependence of the optical depth 
(r) on wavelength. A very rough analytical approximation to 
the extinction law can be written in the form r oc A~", with 
ra ~ 1 — 2 (Draine & Lee 1984). The extinction properties of 
dust particles are still poorly known, both empirically or via 
modelling. As a further complication, it is known that the 
dust obscuration depends on how stars and dust are mixed 
together, so that the extinction law can differ not only from 
galaxy type to galaxy type but even from galaxy to galaxy. 
Besides, little is known about the relative importance of ab- 
sorption to scattering, as well as on the evolution of dust 
in galaxies. Dust is formed during the post-main sequence 
evolution of massive stars, is therefore expected to appear 
in the galaxies soon after the formation of the first stars. 

In the number counts predictions presented here we 
have used a very simple model for dust, which is the same 
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for all galaxy types and at any time (i.e. we do not model 
dust evolution). The extinction law is taken to be oc , 
and the optical depth of an i, galaxy is assumed to be 
r, = 0.3 in the B-band (which, in a standard model, corre- 
sponds to an extinction of Ab ~ 0.15 mag at 2 = 0). We 
also assume that the optical depth depends on the galaxy 
luminosity as r oc LI'^q , to account for the fact that low lu- 
minosity nearby galaxies seem to be less obscured by dust. 
In any case it is worth mentioning that this assumption does 
not have strong effect on the predictions (Campos & Shanks 
1995). As we will comment later, the counts are dominated 
by Lt galaxies, except at the faint limits where the presence 
of the intrinsically faint objects becomes very important. 

The importance of dust in the number counts models is 
mainly due to the combination of spectral red-shifting due to 
the Doppler effect and the shape of the extinction law. The 
effect of dust is more important as the redshift of the galaxy 
increases, because we observe light emitted at shorter wave- 
lengths where the extinction is larger. Thus, the extinction 
by dust has the effect of dimming the brightening of the 
luminosity at high redshifts. In fact, luminosity evolution 
models in which dust is not included overpredict the counts 
at faint limits. 

As can be seen in Figure 3, the go = 0.05 model provides 
a very good fit to the data in ~ 10 magnitudes interval 
range (5 ~ 17 — 27). The go = 0.5 model fails to fit the 
counts fainter than B ~ 24, as it predicts fewer galaxies 
than observed. With respect to the K band, the go = 0.05 
model slightly overpredicts the counts at both the bright 
and the faint end. The go = 0.5 model gives better fit to the 
data down to K ~ 21, but underpredicts the counts fainter 
than that. The discrepancies between models and data in 
the K-band are not very large, and so we regard them as 
due to uncertainties in the spectrophotometric models when 
computing the B — K colour of the model galaxies. 



4 THE REDSHIFT DISTRIBUTIONS OF 
FAINT FIELD GALAXIES 

In order to constrain the number counts models we can make 
use of the redshift distribution N{z) of faint field galaxies as 
an independent test. To compute the N{z) predictions for 
galaxies selected in the magnitude range [mi, 77*2] we just 
have to integrate: 

N(z)dz= / d^A(m,z) (6) 

•J m-i 

In Figures 4 and 5 we plot the N (z) distributions of 
galaxies selected in the B-band to have an apparent magni- 
tude of 21 <B< 22.5 (data taken from CoUess et al. 1990, 
1993) and 23 <B< 24 (data taken from Glazebrook et al. 
1995 and Campos et al. 1995) respectively. In the same Fig- 
ures we also plot predictions from the luminosity evolution 
models (with dust) presented in the previous section. As can 
be seen, the models fit fairly well the data for the shallow 
sample, although for the deep one they predict a high red- 
shift tail which is not observed. The N (z) distribution is 
quite close to the predictions from no-evolution models, also 
shown in the Figures. In fact this came as a surprise, as the 
no-evolution models are unable to fit the counts. 




0.5 

z 



Figure 4. The redshift distribution of galaxies of 21 < B < 
22.5 and predictions from no-evolution and luminosity evolution 
models with dust (LE-j-D). 
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Figure 5. The same as Figure 4 but for 23 < B < 24 galaxies. 
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The discrepancies between data and predictions from 
luminosity evolution models might be due to an observa- 
tional bias rather than a failure of the models (Campos et 
al. 1995). The deep sample is incomplete, i.e. the redshift of 
~ 40% of the objects could not be identified as no spectral 
features were detected. If these galaxies turned out to be at 
high— 2 then the models would fit the observations rather 
well. There is a very simple reason which allows us to sug- 
gest that this might actually be the case. If a galaxy was at 
a redshift higher than 2 ~ 1, then the (quite often strong) 
[011]A3727 emission line would be red-shifted outside the 
optical spectral window (usually running from ~ 4000 to 
8000 A). At shorter wavelengths there are of course other 
spectral features, such as C, Fe or Mg lines, but these are 
usually very weak and so difficult to detect in these faint 
galaxies. We therefore expect that there is a bias against 
the redshift-identification of high-2 objects. A further indi- 
cation pointing out that this might actually be the case is 
found in the colour distribution of the objects. The models 
predict that this evolved high-2 population is bluer, and, as 
shown by Campos et al. (1995), the unidentified galaxies are 
bluer on average than the galaxies with measured redshifts. 

The importance of dust to model the counts is further 
evidenced in the N (z) predictions. The same models pre- 
sented here but without dust (Campos & Shanks 1995) pre- 
dict a much larger redshift tail, totally incompatible with 
the data even in the case that the unidentified galaxies were 
found to be at high-2. 




Figure 6. The number counts in the B and K bands, and pre- 
dictions from a go = 0-5 luminosity evolution (with dust) model. 
The faint end slope of late type galaxies is assumed to evolve as 
a. = a.z=o -\- n ^ z with n = 0.5. 



5 COULD WE FIT THE COUNTS IN A Qo = 0.5 
UNIVERSE MODEL? 

We have shown in the previous sections that simple lumi- 
nosity evolution models which incorporate dust are able to 
provide reasonable fits to the counts and the redshift dis- 
tributions, but only when go = 0.05. As shown in Figure 3, 
the go = 0.5 model underpredicts the counts fainter than 
B ~ 24. 

In general, at any apparent magnitude range the counts 
are dominated by i, galaxies. The contribution of galaxies 
brighter than i,, even if they are placed at higher-z and 
so we survey a larger volume for the same apparent magni- 
tude, is smaller simply because they are less abundant. The 
number density of galaxies fainter than i, is larger, but for 
the same apparent magnitude range they are closer and so 
in a smaller volume. This situation changes as we move to 
lower apparent magnitudes. As clearly illustrated by Driver 
et al. (1994), the bright galaxies start to notice the cosmolog- 
ical volume effect, i.e. the volume decreases, especially in a 
closed geometry, whereas the contribution from faint galax- 
ies situated at lower-2 still increases with a Euclidean slope. 
Therefore at faint apparent magnitudes the counts start to 
be dominated by the intrinsically faint galaxies. 

In order to reconcile the predictions from the go = 0.5 
model with the observations we could introduce a population 
of intrinsically faint objects (eg. dwarf galaxies, low surface 
brightness galaxies...) in the LF. In fact we know that dwarfs 
are very common in our local neighbourhood, e.g. in local 
clusters of galaxies the percentage of dwarfs is as large as 
~ 30 — 50% (Binggeli, Sandage & Tamman 1985; Phillips et 
al. 1987; Davies et al. 1988). 



If we steepen the faint end slope of the LF in order to 
fit the counts with the go = 0.5 model then we would also 
predict a large population of low-z objects in the N{z) dis- 
tributions which is not observed. A possibility is to assume 
that dwarfs were more numerous in the past, i.e. their num- 
ber density increases toward higher redshifts. In Figure 6 we 
plot the number counts in the B and K bands together with 
predictions from the go = 0.5 luminosity evolution model 
(with dust). The difference between this model and that 
presented in section 3 is that here we assume that the faint 
end slope of the late type galaxies LF evolves with the red- 
shift as a = az=o + n * z with n = 0.5. The predictions are 
now much closer to the observations. On the other hand, 
the N{z) predictions of this model are quite similar to those 
shown in section 4, as the contribution of dwarfs starts to 
be noticeable at the very faint end, beyond the limits of 
the deep redshift surveys. Nevertheless it is worth mention- 
ing that the models shown in Figure 5 for the 23 <B< 24 
sample slightly overpredict the number of galaxies at low-z 
(z ~ 0.1 — 0.3). We suggest that a small percentage of the 
unidentified galaxies might be low-z dwarfs, where the star 
formation has recently been stopped and therefore they are 
still blue but showing featureless spectra (Campos 1995). 

The evolution of the faint end slope of the late type 
galaxies is a rather ad hoc ingredient which is based on the 
models to explain the nature of dwarf galaxies. Dwarf galax- 
ies are thought to form stars in intermittent bursts. After a 
burst of star formation is switched on, the small mass of the 
galaxies might not be able to retain the gas supply against 
the formation of the galactic wind that might follow the 
overlap of supernova explosions (Dekel & Silk 1986). If the 
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Figure 7. The I-band counts of Glazebrook et al. (1995) from the 
HST Medium Deep Survey, for ellipticals plus spirals (triangles) 
and irregulars (circles). We also show predictions from the go = 
0.5 model with evolving faint-end slope. 

star formation is abruptly stopped, then the luminosity of 
the small galaxies will fade away. Therefore if we include the 
population of dwarfs in the late-type group of galaxies we 
must allow for luminosity-dependent luminosity evolution, 
which we have tried to mimic by changing the slope of the 
LF. 

Using the Hubble Space Telescope HST we are now able 
to obtain images of very high resolution, such that we can 
resolve the morphology of the faint field galaxies. There- 
fore it is possible to obtain the number counts for galax- 
ies divided according to their morphology. Recent work on 
number counts based on HST images (Glazebrook et al., 
1995; Driver et al. 1995) has shown that the observed late- 
type/lrregular population is as much as a factor of ~ 10 in 
excess of the conventional no-evolution models. This excess 
has been interpreted by these authors as due to a popula- 
tion of dwarfs that has either faded away, or has merged. 
In Figure 7 we plot the HST counts for early-intermediate 
galaxies and late/Irregulars (data taken from Glazebrook et 
al.) together with predictions from the go = 0.5 model with 
evolving LF slope. The model fits the data quite well, ex- 
cept for the last bin where it overpredicts the number of 
early /intermediate galaxies while underpredicting the num- 
ber of Irregulars. We regard this discrepancy as possible 
misidentification of the high-2 progenitors of present day 
ellipticals/spirals or just to the simplicity of the model. 

6 CONCLUSION 

Number counts is a conceptually simple method of measur- 
ing the value of go . However before we can properly address 



this cosmological problem we have to solve an astrophysi- 
cal problem concerning the formation and evolution of the 
galaxy population. This is proving to be a difficult task, 
as many different processes, such as evolution of the galax- 
ies luminosity or mergers (see e.g. Guiderdoni & Rocca- 
Volmerange 1991) are involved. 

Even if the modelling might seem a rather complicated 
problem, we have shown that simple models that consider 
the evolution of the galaxies luminosity and the extinction 
produced by dust are able to provide reasonable fits to both 
the number counts and the redshift distributions of faint 
field galaxies. 

These models only provide a fit to the counts at the 
faint end when go = 0.05. For go = 0.5 the models un- 
derpredict the counts. A possibility of reconciling models 
and observations is by means of a population of intrinsically 
faint objects which becomes increasingly more important at 
higher redshifts. 
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